A0AO2176O 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


[s IMULTANEOUS  IDENTIFICATION  OF  SHORT  PERIOD  AND  PHUGOII 
f PARAMETERS  USING  AN  ADVANCED  LIKELIHOOD  METHOD 


Frederick  T.  Bryan 
December  1975 


Thesis  Advisor  Ronald  A.  Hess 

Approved  for  public  release;  distribution  unlimited. 


UNCLASSIFIED 

security  cl »»unc nio*  or  tmh  pace  (v*««i  o««  tmm,m€> 


REPORT  DOCUMENTATION  PAGE 


I.  REAORT  NUMIIr 


i*.  OOVT  ACCESSION  HO 


TlTl  r f"J 


SIMULTANEOUS  IDENTIFICATION  OF  £HORT  PERIOD  AND] 
JHUGOID  STABILITY  PARAMETERS  USING  AN^DVANCED 
XlKELIHOOD  METHOD, 


£ 


Frederick  T.^ryanjt  Major,  USMC 


4.  ecnronuiNO  oroanization  name  ano  adorcii 


Naval  Postgraduate  School 
Monterey,  California  93940 


ii.  controlling  ornce  name  and  address 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


».  reorient'*  catalog  numoca 


t.  tyre  or  report  • re  nop  covered 


Aeronautical  Engine^r- 
December  1975 


*•  AERrOAMINO  OAO.  REPORT  NUMBER 


&Sr  5 4-UeS\S 


10.  AROQRAM  ELEMENT.  PROJECT,  TASK 

* 1IT  f 


AREA  « VORK  UNIT  NUMBERS 


Naval  Postgraduate  School 
Monterey,  California  93940 


14  moniTOAinO  AGCnCV  nauI  • AOORESSdJ  dlllmrmnt  from  Controlling  OHIcm) 


Naval  Postgraduate  School 
Monterey,  California  93940 


120 


l».  SECURITY  CL  A**,  (ml  thlm  tipmrt) 


Unclassified 

lassiai 

schedule 


I la.  OECL  ASSlAl  CATION/ OOVNO  A ADlMO 

ED  ■■  - 


l«.  DISTRIBUTION  STATEMENT  (ml  inlm  Kmpm/t J 


Approved  for  public  release;  distribution  unlimited 


17.  Distribution  »T  ATCmEmT  (ml  ihm  mmnrmcl  mnlmrmg  In  BlmmO  70,  II  Olllmrmnt  Irmm  Kmpmrt) 


l».  SUAALEMEnT  A AY  notes 


IS.  KEY  VO  AOS  (Cool  Invm  on  ,m  mmrmm  mlOm  II  nmcmmmmry  an  a itfooii  ty  Or  Wo  as  nivloj 


Maximum  Likelihood  Identification 
Stability  Parameters 
Stability  Derivatives 
Flight  Test  Data  Reduction 


10 


ACT  fCmiimif  on  It  m4  ftfmffff  kf  lltll  mmtm) 


An  investigation  was  conducted  to  determine  the  feasibility  of  obtaining 
short  period  and  phugoid  stability  derivatives  from  one  maneuver, 
simultaneously  . It  was  concluded  that  the  uiaxiraum  likelihood  identification 
program  SCIDNT  I showed  great  promise  in  obtaining  the  short  period  and 
phugoid  stability  derivatives  from  one  maneuver.  Extraction  of  the  short 
period  stability  parameters  in  the  presence  of  the  phugoid  was  easy,  straight 
forward,  and  yielded  results  similar  to  those  obtained  from  pure  short  period- 


‘lift 


0D  1473  A EDITION  or  I NOV  «•  II  OBSOLETE 


Pd 


(Page  1) 


S/N  0 10  2*  0 14*  410  1 | 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  or  THIS  RAO!  (•*>•"  SMI  On  I 

&r/¥s'V  f 


A 


UNCLASSIFIED 


Jucuwity  cl aSSiFicatiom  or  this  n»i« 


•tUV 


data.  Estimation  of  the  phugoid  stability  parameters  At&s  possible  when  they 
were  estimated  in  conjunction  with  Z^,  hU,  and  €L  . It  was  recorTimended 
that  a new  set  of  data  be  obtained  at  a flight  condition  where  the  phugoid  is 
at  least  moderately  damped  and  that  this  data  be  analyzed  to  resolve  the 
present  anomalies. 


w 


o' 


vs>*- 


_ 1 


DD  FornQ  1473 
S/rJ  01&2-014-G601 


**CU«ITV  CU*»*l'IC*TlOK  of  ThI»  Du* 


SIMULTAN IOUS  IDENTIFICATION  OF  SHORT  PERIOD  AND  PHUGOID 
STABILITY  PARAMETERS  USING  AN  ADVANCED  MAXIMUM  LIKELIHOOD 

METHOD 


by 


Frederick  T.  Bryan 
Major  United  States  Marine  Corps 
E.S.B.A,.  Boston  College  1963 
M.S.A.E.,  Naval  Postgraduate  School,  1974 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


AERONAUTICAL  ENGINEER 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  1975 


Author: 




Approved  ty  = <R Q. 


Hies  Is  X3  visor 


3 


ABSTRACT 


^ ■ 


L 

F 

f. 

I 

\ 


li 

r 

I 


An  investigation  was  conducted  to  determine  the 

feasibility  of  obtaining  the  short  period  and  phugcid 

stability  derivatives  from  one  maneuver, 

“simultaneously".  It  was  concluded  that  the  maximum 

likelihood  identification  program  SCIDNT  I showed 

great  promise  in  obtaining  the  short  period  and 

phugoid  stability  derivatives  from  one 

maneuver.  Extraction  of  the  short  period  stability 

parameters  in  the  presence  of  the  phugoid  was 

easy,  straight  forward,  and  yielded  results  similar  to 

those  obtained  from  pure  short  period  data.  Estimation 

of  the  phugoid  stability  parameters  was  possible  when 

they  were  estimated  in  conjunction  with 

Z , X , M , and  0 . It  was  recommended  that  a new  set 
o o o o 

of  data  be  obtained  at  a flight  condition  where  the 
phugoid  is  at  least  moderately  damped  and  that  this 
data  be  analyzed  to  resoivo  the  present  anomalies. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Plight  testing  of  military  aircraft  has  always  been  a 
time  consuming  and  expensive  proposition.  The  relatively 
recent  introduction  of  the  new  Flying  Qualities 
Specification  (MTL-F-6765B)  has  introduced  many  new 
parametric  reguir ements  which  reguire  that  many  of  the 
aircraft  airframe  parameters  be  determined  to  insure 
specification  compliance.  As  a result  of  these  new 
requirements  the  Naval  Air  Test  Center  (NATC)  initiated  a 
research  program  to  develop  airframe  parameter 
identification  technology  for  use  in  flight  testing  Davy 
aircraft.  This  program  was  initiated  in  1971  between  the 
Naval  Air  Systems  Command  (NAVAIR)  and  NATC.  At  the  same 
time  the  Office  of  Naval  Research  (ONR)  issued  a contract  to 
Systems  Control,  Inc.  (SCI)  to  advance  the  state-of-the-art 
in  parameter  identification.  In  1973  a concentrated  effort 
was  begun  by  the  Navy  to  continue  the  development  of  this 
new  technology  by  forming  a joint  program  between 
NATC,  NAVAIR,  ONR  and  SCI.  As  a result  of  this  development  a 
digital  FORTRAN  computer  program,  SCIDNT-I  (hereafter 
refered  to  as  SCIDNT)  was  written  for  NATC  to  extract 
aircraft  stability  and  control  coefficients  from  flight  test 
data.  Ref.  1.  It  was  hoped  that  by  using  this  computer 
program  it  would  be  possible  to  significantly  reduce  the 
flight  test  time  in  the  air,  and  the  ground  data  reduction. 
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B.  PURPOSE 


The  purpose  of  this  research  was  to  develop  parameter 
identification  test  methods  for  determining  aircraft  phugoid 
and  short  period  characteristics  from  one  aircraft 
maneuver.  This  research  further  entailed  the  determination 
of  proper  aircraft  maneuvers/pilot  inputs  and  detailed  data 
analysis  proceedures  using  the  advanced  maximum  likelihood 
identification  technique. 


C.  DESCRIPTION  OF  AIRPLANE 

The  test  airplane  was  demonstration  F-14A  ship  number  8 
BU.  NO.  157S87.  The  airplane  was  essentially  typical  of  a 
production  F-14A  airplane  relative  to  these  tests.  An 
aircraft  description  is  contained  in  Appendix  A.  photographs 
and  a three  view  drawing  of  the  test  airplane  are  shown  in 
Appendix  A , Figures  1,  2 and  8. 

D.  MAXIMUM  1IKELI HOOD  PARAMETER  IDENTIFICATION  THEORY 

A widely  used  rule  for  choosing  the  sample  statistic  to 
be  used  in  parameter  estimation  is  called  the  method  of 
maximum  likelihood  Ref.  2.  This  method,  in  addition  to 
prescribing  the  statistic  which  should  be  used,  provides  an 
approximation  of  the  distribution  of  that  statistic  so  that 
approximate  confidence  intervals  may  be  constructed.  The 
method  of  maximum  likelihood  makes  use  of  the 
sample- likelihood  function.  Assume  we  observe  or  record  a 
set  of  observations  on  a system,  whose  model  has  p unknown 
parameters  6.  -We  then  can  write  the  joint  probability 
distribution  of  a random  sample  for  any  given  set  of  values 
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of  the  parameters  6 from  a feasible  set  8.  He  then  can 
assign  a probability  p(z|0)  to  each  outcome  z.  If  ve 
observe  or  measure  the  outcome  of  an  actual  experiment  to  be 
z we  then  mould  like  to  know  which  set  or  sets  of  values  of 
6 might  have  resulted  in  these  observations.  This  concept  is 
contained  in  the  likelihood  function  aC  (9|z).  This  function 
is  of  fundamental  importance  in  estimation  theory  because  of 
the  likelihood  principle  of  Fisher  and  others  which  states 
that  if  the  system  model  is  correct,  all  information  about 
unknown  parameters  is  contained  in  the  likelihood 
function.  The  maximum  likelihood  estimator  of  8 is  the  value 
8 which  causes  the  likelihood  function  <£(61  z)  to  be  a 
maximum  [ 2 J . 


d = max  (@  I z) 
8C« 


From  a purely  mechanical  viewpoint  this  method  can  be 
' . summarized  as  follows: 


Find  the  probability  density  functions  of  the 
observations  or  measured  values  for  all  possible 
combinations  of  the  unknown  parameter 
values.  Select  the  density  function  whose  value  is 
highest  among  all  density  functions  at  the 
measured  values  of  the  observations.  The 
corresponding  parameter  values  are  the  maximum 
likelihood  estimates. 


Assume  8 can  take  three  possible  values:  8^,  6 , and 

6^.  Further  assume  that  the  probability  density  functions  of 

observations  z for  these  three  values  of  8 are  as  shown  in 
Figure  1.  Then,  if  the  actual  observation  is  z,  8 is  the 

maximum  likelihood'  estimate  of  0. 
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Figure  1. 


z 

observed 

Maximum  Likelihood  Estimates 

Owing  to  the  monotonic,  one-to-one  relationship  between 
the  likelihood  function  and  its  logarithm,  i.e.  they  have  a 
maximum  at  the  same  value  6,  it  is  more  convient  tc  work 
with  the  negative  of  the  logarithm  of  the  likelihood 
function.  Starting  from  a priori  values,  the  parameters  are 
updated  so  that  the  value  of  the  negative  log- likelihood 
function  J at  the  observed  value  of  the  outputs  decreases 
continuously. 

As  sample  statistics  or  functions  of  the 
observations,  the  estimators  are,  before  the  experiment, 
random  variables  and  can  be  studied  as  such.  In  fact,  the 
major  advantage  of  maximum  likelihood  estimators  is  that 
their  properties  have  been  thoroughly  studied  and  are  well 
known,  virtually  independently  of  the  particular  density 
function  under  consideration.  Additionally,  the  maximum 
likelihood  method  can  be  used  with  linear  or  nonlinear 
models  in  the  presence  of  process  and  measurement  noise.  The 
estimator  of  0 is  approximately  normally  distributed  for 
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large  n,  and  is  asymptotically  unbiased,  efficient, 
consistent,  and  sufficient,  [21. 

The  details  of  the  maximum  likelihood  identification 
proceedure  are  presented  in  the  following  section. 

! 

i 

4 

E.  MAXIMUM  LIKELIHOOD  IDENTIFICATION  OF  PARAMETERS 

t 

■ The  derivation  follows  closely  that  presented  in  Ref. 

J 3,  and  is  presented  herfein  for  completeness  and  reader 

convenience.  The  initial  eguations  will  be  concerned  with 
general  nonlinear  systems. 

Consider  the  general  nonlinear  aircraft  equations  of 
> motion 


i = £(x,  u,  e,  t)  ♦ r«>.  t) 


0 < t < T 


( 1 ) 


E(x(0))  = X (6) 

o 


E{(x(0)  - x (6))  (x  (0)  - x (9))  } = P (9) 


where 


( 2 ) 


£ 


- 

i 

P\ 


X(t) 

is 

n 

X 

1 

state  vector 

u<t) 

is 

1 

X 

1 

input  vector 

e 

is 

P 

X 

1 

vector  of  unknown  parameters 

P(e,t) 

is 

n 

X 

9 

process  noise  distribution  matrix 

w(t) 

is 

9 

X 

1 

random  process  noise  vector 

Sets  of  m measurements  y (t  ) are  taken  at  discrete  times  t 

k k 
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( 3 ) 


v (t)  and  v (t  ) are  Gaussian  random  noises  [2],  with  the 
Ic 

following  properties 


E(w(t))  = 0 E (v  (t  ) ) = 0 

k 

T 

E(w(t)  V (t  ))  = 0 
k 

T 

E(w(t)  W (T))  = Q (6,  t)  6 (t-T  ) 

T 

2 (v  (t  ) v (t  ) ) = R (9,t  ) <?  ( U ) 

j k j jk 

The  unknown  parameters  are  assumed  to  occur  in  the 

functions  f and  h and  in  matrices  P,  Q,  fi,  P , and  x . In 

1 o o 

the  following  analysis,  the  model  and  the  functional  form  of 

f and  h is  assumed  known  correctly. 

The  set  of  observations  y (t  ),y(t  y (t  ) 

» a n 

constitutes  the  outcome  Z in  this  case.  The  likelihood 

function  for  B,  which  has  the  same  form  as  the  probability 
of  the  outcome  z for  a certain  value  of  parameters  6,  is 
given  by 

^ (6|z)  = p(z<6) 

fc  P(y  <t  ) , y<t  ) ,...,y(t  ) |6) 

• ^ N 
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< 5 ) 


*(©|z)  = P (Y  16) 
N 

where 


T « (y(t  ) ,...,y  (t  )} , k = 1,  2,  • • • # H 
k 1 k 

p(Y  |0)  = p|y  <t  ) < Y ,©)•  p(Y  |6) 

N N N®1  N® 1 

* '"'VVl'  61  P(y(t»-l"I»-2'  *>  «V!ie| 

which  by  successive  application  of  Bayes  rule,  [2],  becomes 

N 

p{Y  |0)  = TT  P(y  (t.)  |Y.  , 0)  ( 6 ) 

N 1 1 11-1 

i=1 


The  log- likelihood  function  is 


log  <^<6|z))  = Y, 


1-1 


log  {p  (y  (t.)  |Y.  , 

i i-I 


0) ) + constant  ( 7 ) 


To  find  the  probability  distribution  of  y(t  ) given  Y and 

i i-1 

0,  the  mean  value  and  covariance  are  determined  first. 


E(y<V  1 0)  & y (i I i-1) 


( 8 ) 


The  expected  value  or  the  mean  is  the  best  possible  estimate 
of  measurements  at  a point  given  the  measurements  up  to  and 
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including  the  previous  point. 


cov(y(t.)|T  ,0)  = B{(y(t  ) - J (i  I i- 1) ) (y  (t .)  - / (i|  i-T) ) 
a.  i-1  i i 

V"(i)T} 

6B(i)  ( 9 ) 


V'  <i)  are  called  the  innovations  at  point  i and  B(i)  is  the 
innovations  covariance.  Since 

y (t.)  - E(y  (t.)  i y . ,e)  = \^{i)  ( io  ) 

i ii-i 

it  follows  that  y<t  } given  Y and  0 have  the  same 

i i-1 

distribution  as  (i) . It  has  been  shown  that  as  the 

sampling  rate  is  increased,  the  innovations  V (i)  tend 
towards  having  a Gaussian  density.  Assuming  a sufficiently 
high  sampling  rate,  the  distribution  of  V'  (i)  and, 
therefore,  y (t  ) given  Y and  6 is  Gaussian,  i.e.. 


p(y(V  ,Vi'e) 


1 

exp 

2 

■ 


V'(i)  B-MD  v'  (i) 


, b/2  i 1/2 

(2TT)  I B (i)  | 


( 11  ) 
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log  {p(y(tA)  / <d  B-i(i)  ^(i) 


log|B(i)j  ♦ constant 


( 12  ) 


The  log-likelihood  function  of  equation  ( 7 ) can  be  written 


H 

log  ( *(612)}  = - — £ 


l*  (i)  b~  1 (i)  v"  (i) 


♦ log  |B  (i)  |} 


( 13  ) 


An  estimate  of  the  unknown  parameters  is  obtained  by 
maximizing  the  likelihood  function  or  the  log- likelihood 
function  from  the  feasible  set  of  parameter  values. 


6 = ef«  109  ( 


( 14  ) 


= max 


r i ” 

[ ~ ~ X (i)  (i)  / (i)  * log  I B ( i)  j}  1 

i=  1 J 


( 15  ) 


The  log-likelihood  function  depends  on  the  innovations 
and  their  covariance.  To  optimize  the  likelihood  function,  a 
way  must  te  found  for  determining  these  quantities.  Both 
innovations  and  their  covariance  are  outputs  of  an  extended 
Kalman  filter. 


The  extended  Kalman  filter  is  conventionally  divided 
into  two  parts.  In  the  first  part,  called  the  prediction 
equations,  the  state  equations  and  state  estimate  covariance 
equations  are  propagated  in  time  from  one  measurement  point 
to  the  next.  In  the  second  part,  called  the  measurement 
update  equations,  the  measurements  and  associated 
measurement  noise  covariances  are  used  to  improve  state 
estimates.  The  covariance  matrix  is  also  updated  at  this 
point  to  reflect  the  additional  information  obtained  from 
the  measurements. 

1 . predict  ion  Equations 

The  state  prediction  is  done  using  the  equations  of 

motion.  Appendix  B.  Starting  at  time  t with  current 

i~1 

estimate  2(i-1|i-1)  of  the  state  x (t  ) and  the  covariance 

i 

P(i-1|i-11,  the  following  equations  are  used  to  find  the 

predicted  state  x(i|i-1)  and  the  associated  covariance 
P (i | i-  1)  . 


7-  x(t|t.  ) = f(x(t|t.  ),u(t),e,t)  ( 16  ) 

dt  i-l  i-i 


P(t|t.  ) = F(t)  P(t|t.  ) ♦ P(t|t.  ) P (t) 
1-1  1-1  1-1 


♦ r^f1  < 17  > 


t s t < t 
i-1  i 


2a 


The  n x n matrix  P is  obtained  by  linearizing  f 
about  the  best  current  estimate 


F(t) 


jf(x(tit.  ) ,u(t)  ,e,t) 

i- 1 


Dsing  { 16  )-(  18  ),  we  can  obtain 


x(t.|t,  ) £>  x(i|i-1) 

i i-l 


( 18  ) 


and 


P(t.  it.  ) 6 P(i|i-D 

i i-l 


( 19  ) 


Thereafter,  the  measurement  update  eguations  are  used. 

2*  Measurement  Update  Equations 

The  covariance  and  state  estimate  are  updated  using 
the  measurements.  The  necessary  relations  are  derived  in 
Ref.  4,  and  are  presented  here  without  proof.  The  innovation 
and  its  covariance  are 


Y^i)  = y (i)  - h(x(i|i-1),u(t  ),e,t.)  ( 20  ) 

i i 


B (i)  = H (i)  P (i|  i-1) 


T 

H <i> 


♦ R 


( 21  ) 
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where  H is  obtained  by  linearizing  h. 


H(i) 


jh{x(i|i-1)  ,u(t  ,©,t  ) 
1 i 

ax  <i|  i- 1) 


( 22  ) 


The  Kalman  gain  and  the  state  update  equations  are 


T 

K(i)  = P(ili'l)  a {i)  B-»(i)  ( 23  ) 

x ( i 1 i)  = x (i  | i—  1 ) ♦ K (i)  v'(i)  ( 24  ) 


Finally,  P{iji),  the  covariance  of  error  in  updated 
state,  is  obtained  by 


P(i|i)  = (I-K(i)  H (i) ) P (i  | i-1) 


( 25  ) 


One  can  now  return  to  the  time  update  or  prediction 
equations. 

3.  Optimisation  Procedure 

Many  possible  numerical  procedures  can  be  used  for 
this  optimization  problem.  Modified  Newton-Raphson  or 
Quasilinearization  hare  been  found  by  experience  to  give 
quicker  convergence  than  most  procedures  such  as  the 
conjugate  gradient  or  the  Da/ison  method.  The  modified 
Newton-Raphson  method  is  a second  order  gradient  procedure 
requiring  computation  of  first  and  second  order  parti als  of 
the  log-likelihood  function. 
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also 


Hi*  T 

a*  log  < <(©U))  y r jV 

^ ‘ ft  L 39 


<iLB-Mi)  av'(i) 


a© 


a v''  (i) 


a© 


j 


fle 


j 


8B(i) 

B-I  (i)  — — B-l  (i)  V'(i) 

a®. 


^ t . aB(D  ae(i)  , . 

*v  (i)  B” 1 (i)  — B->(i)- B~*U)V'(i) 


- — Tr  J B-J  (i) 


a« 


dB  (i) 


a© 


a© 


B-»  (i) 


j 


k 

jBgj 

*** 


+ B-i  (i)  V' (i) 


- V'T(i)B-*{i)  B-»  (i)  V (i) 


* — Tr  B-»  (i) 
2 1 


a23g)  V 

3W. 


( 27  ) 


j ~ 1#2»  « • t »p 


The  last  three  terms  in  the  equation  for  the  second 
partial  of  the  log- likelihood  function  involve  second 
partials  of  innovation  and  its  covariance.  These  terms  are 
usually  dropped  which  yields  the  second  partial 
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approximation 


a*  log  (se(9|z)) 


Nr  ^ 

•sh= 


(i)  . . a v'(i) 

1 — B~ 1 (i) 


a® 


j 


a/  (i)  jB(i)  , 

-B-Mi)-— B-MDV'U) 


*k 

a ^(D 


a® 


j 


a® 


j 


as  (i>  ^ . 

B-i  (l)  — B~ 1 (1)  V (1) 

a®. 


.#  t as(i)  ae(i) 


Tr  B-  » (i) 


a® 


as(j) 

a®, 


a® 


B~ 1 (i) 


k 

aB(i) 

a® 


k /j 


( 28  ) 


From  equations  20  and  21  we  can  obtain  the  gradients 
of  innovation  and  its  covarience  for  parameter  6 as# 


a /U) 


a® 


j 


x=x (i | i'  1) 


ax(iU-P 

a® 


a_h 

a® 


J = 1#2#...#p 
i * 1, 2 , . . . # N 


( 29  ) 
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i 


iiiiL  = iiiiiL  p(i| i-i>  HT(i»  . ■ 

J°i  J8J  J9l 


an  (i)  jr 
+ Hli)  PUIi-1>— — 

j j 


( 30  ) 


j * • /P 

i = 1,2,.  . . , H 


Recursive  equations  can  be  obtained  for  gradients  of 
the  predicted  state  and  its  covariance.  This  is  done  iD 
stages  by  using  the  prediction  and  measurement  update 
equations  of  the  Kalman  filter.  Differentiating  equations 
16-18  with  respect  to  0 yields 


Jx(tJt.  ) Jf  <x  (t  1 1 . ) , u (t)  , t,0) 

d 1-1  i-l 


dt  J0 


Jf(x(t|t,  ),u(t),t,0)  Jx(t|t  ) 

i-l  i-1 

’ a2(tit."  > * &T 

i-i  j 


. ax  (9) 

ax  (oio)  o 


( 31  ) 


. JP(tjt.  ) 
d 1-1 

Tt  a«  — 


jp(tit  ) 

ap(t)  i-v 

ae  i-i  ae 


ap  tit  ) t 

i-1  T jF 

♦ P (t)  ♦ P(t|t  ) 

ae_j  'i-i  ao^ 


ap  («) 

ap(oio)  o 


t 5 t it 
i-1 


( 32  ) 


a2f  (x(t|t.  ) ,u  (t)  ,0,t) 

3F  (t)  = i^l 

J0 . ~ ae  ax(t|t.  ) 

j j x-1 


( 33  ) 


j = 1,2, 


The  sensitivity  functions,  used  in  the  SCIDNT  integration 
routines,  are  updated  at  aeasurement  points  by 
differentiating  eguations  22  - 25  with  respect  to  8^. 


a*b(x(i|i-l),u(t  ),8,t  ) 

g H ( x } x X 

a« _ 3 ae  asum-D 


< 34  ) 
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aK<i)  ap  (i 1 1- 1)  t 

0 — — H (i)  B*1  (i) 


a© 


j 


a© 


j 


T 

jh  <i) 

♦ P (i|  i-1)  — B~ 1 (i) 


a© 


j 


T JB  (i) 

- P (i|  i-1)  H (i)  B-i  (i)  — - — B-‘  (i)  ( 35  ) 


a© 


j 


ax(i|i)  j)x(i|i-1)  3K(i)  ^ . . av'U)  _ 

♦ /(i)  + K(i)— ( 36  ) 


as 


j 


a© 


j 


a© 


j 


a© 


j 


aP(ili) 

a©. 


= (I-K(i)  H (i) ) 


apaii-i) 

a©. 


aMi> 


a© 


H(i)  P(i|i-1) 


j 


ami)  . . 

- K (1)  — — P(i|l-1) 


a© 


j 


( 37  ) 


j ~ )#2#3/..vrp 


Tho  negative  of  the  matrix  of  second  partials  of  tne 
log-likelihood  function  is  called  the  information  matrix 
M.  The  step  size  A 9 for  parameter  estimates  is  given  by 


....  a log  ( ae(©U) 

A 9 * . 

a© 


( 38  ) 
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The  information  matrix  provides  a lover  bound  on  parameter 
estimate  covariances,  i.e. , 


E (0-d) 


(0-0)  £ H-i 


( 39  ) 


This  is  the  Cramer-Bao  lover  bound.  The  maximum  likelihood 
estimates  approach  this  bound  asymptotically. 

4 . Linear  Systems 

In  a linear  system,  the  functions  f and  h are 
defined  as 


a 

j 

i 

3 

1 


* 

-J 

i 

i 


f(x,u,e,t)  £F(e,t)x  *■  G(9,t)u 


and  h (x,u,9,  t)  £>H(0,t)x  * D(0,t)u 


( 40  ) 


Por  the  linear  case,  it  has  been  shown  that  if  the 
model  is  correct,  the  innovations  are  white  and  have 
Gaussian  density  at  the  true  values  of  the  parameters.  The 
assumption  of  fast  sampling  rate  is,  therefore,  not 
necessary , 

The  basic  algorithm  for  solution  of  the  linear  case, 
is  the  same,  however,  some  of  the  equations  can  now  be 
simplified.  The  equivalence  between  F(t)  and  H(i)  of  { 18  ) 
and  ( 22  ) and  F(0,t)  and  of  ( 40  ) is 

obvious.  Equations  ( 16  } and  ( 20  ) now  become 


It 


f (e,t>  xttit^)  ♦ G(e,t)  u(t) 

{ 41  ) 


and 

v'(i)  = J(i)  - H(9,t.)  x (i|  i-1)  - D(0,t.)  u(t.) 

1 xx 

( 42  ) 


Equations  ( 29  ) and  { 31  ) can  be  written  as 


and 


3 V'(i) 


a« 


j 


3* (i| i-1) 

H (6/  t ) 

1 8® 


3H (0,^) 

— x (i|  i-1) 


a® 


j 


ao<e,t.) 

U(t  ) 

ae3 


( 43  ) 


d ax<t|Vi} 

dt  a®. 


= F (6  / 1) 


axitit.  ) 

i-i 


a® 


j 


aF(®,t;  a _ . 

♦ — *<t|t.  ) 

a®  i-i 


aG<e,t)  _ 

♦ u(t) 


a® 


j 


( 44  ) 


All  other  equations  reiaain  the  same. 

There  is  a considerable  reduction  in  the  computation 
requirement  for  time-invariant  linear  systems.  In  this 
case,  matrices  F,  G,  H,  D,  p , Q and  H and  their  derivatives 
with  respect  to  the  parameters  are  constant. 
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5.  yipe  Invariant  Linear  Systems  In  Statistical  Steady 
State 

In  many  aircraft  applications,  the  Kalman  filter  is 
in  steady  state  for  the  duration  of  the  experiment.  This 
occurs  when  the  Kalman  filter  is  in  operation  for  a 
sufficiently  long  time  and  the  process  and  measurement  noise 
covariances  are  constant.  The  Kalman  gain  and  the 
innovations  and  the  state  covariances  approach  constant 
values.  The  time  update  ana  measurement  update  equations  for 
the  covariances  are 


~p<tit.  ) = pp(t|t.  > ♦ p(tu.  )ft  ♦ rorr 
dt  r-1  1-1  1-1  1 1 

( 45  ) 

T 

K = P (i | i—  1 ) H B~  1 ( 46  ) 

T 

B = HP  {i  | i- 1)  K + R ( 47  ) 

P(iji)  = (I-KH)  P (if  i- 1)  ( 48  ) 

By  definition  of  the  steady  state 
P (i-1|i-1)  = P (i I i) 

if' 


Therefore,  from  ( 45  ) 


T 

P At  P At 

P (i|i-1)  = e P (i-1|i-1)  e 


r 


P(t.-r)  nnT  p (t.-  t) 

j Q I e x 4T 


e x 


i-1 


F At 

= e 

t 


(I-KH)  P (i | i- 1 ) e 


T 

F At 


F(t  -T  ) __T  F (t  -T) 

e i Pop  e i dT: 


t 

i-1 


&0(At)  (I-KH)  P (i  | i- 1 ) *S(  a t)  «•  Q« 

T 

= 0{At)  (P  (i  | j-1)  - KBK  ) tf(At)  + Q' 

( «9  ) 


Using  ( 46  ) , ( 47  ) and  ( 49  ) , we  can  solve  for  P(i|i-1) 
and  then  find  K and  B.  Also,  it  can  be  shown  that 


iP 


JL 

1 je 


T 

A + 
1 


A 

2 


T 

- ftPA^Pj# 


( 50  ) 


tlK  3P  T T 

* (I-KH)  H (HPH  ♦ B)  “ * + A ( 51  ) 

a6j  2 
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r 


it  jjp  t an 

— PH  ♦ H -JL—  H + HP ♦ 

i as  ae  3® 

j 3 3 3 


( 52  ) 


A = f!  (I-KH) 


w t a* 

& = — — (I-KHJP0  ♦*Z5(I-KH)P  

2 ae  a® 


a«  t ao' 

- 0K P0  * — 

ae  a» 

j 3 


aH  t aa  t 

A = B~ 1 H ♦ H 3“ 1 PH 

3 ao  3d 

3 3 


aH  aa 

♦ HP  — ♦ - — b-»H 

ae  ae 

3 3 


( 53  ) 


p & P(ili-I) 


Thus,  it  is  possible  to  solve  for  using  ( 50  ) and  then 

JK  fltB 

find and  from  ( 51  ) and  ( 52  ) . Equation  ( 50  ) 

j»  je 

jP 

is  a linear  equation  in  • and  the  coefficient  of  the 

unknown  matrix  does  not  depend  on  the  parameter 


0^.  Thus,  the  sensitivity  of  the  state  covariance  matrix  can 

be  determined  very  quickly  for  all  parameters.  Once  the 
sensitivity  of  P,  K and  B for  unknown  parameters  is 
determined,  only  state  sensitivity  equations  need  to  be 
updated.  The  computation  of  state  sensitivity  functions  can 
be  reduced  to  many  fever  equations. 

An  approximation  simplifies  the  problem  further.  The 
unknown  parameters  are  defined  to  include  elements  in  K and 
B matrices  instead  of  Q and  R.  Optimizing  the  log- likelihood 
function  for  parameters  in  B gives 


N 

« 1 r*»  , . .T  . 

B = — ) 'Z  (i)  / (i) 

i=  1 


( 54  ) 


The  gradient  of  the  log- likelihood  function  with  respect  to 
other  unknown  parameters  is 


<]  Log  ( of  (6 1 z)  } 


A 


a v*  (j) 

a® 


( 55  ) 


The  sensitivity  of  innovations  to  parameters  is  determined 
using  the  following  recursive  equations 


d . 

— x(t|t.  ) = Px(t|t.  ) ♦ Gu  (t) 
dt  1-1  i-l 
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1 i 

■ f 


d i-1  JP  A 

= _S x (t  1 1 ) ♦ F 

dt  a®  i-i  a® 


j 


j 

♦ u(t) 


j 


j — If 2f axfP  ^ S t S t 


( 56  ) 


v'(i)  = y (i)  - Hx(iii-I)  - Du<t.) 

i 


x (i  j i)  = x <i  1 i—  1 ) + K V"(i) 


( 57  ) 


a v'U) 


a® 


an  A . . aMi|i-i)  3D 

x (i  | i—  1 ) - H - ; u(t.) 
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a*  (i|i) 
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j *1,2, 


ax(ili-l)  3K  . . a V'(i) 

* V'(i)  + K 


a® 


a® 


j 
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i. 
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Note  that 


il 

a® 


= 0 if  6 is  not  an  element  of  K matrix 

j 


( 59  ) 


= i if  e & k 

j • k • j j'k' 
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where  I is  a matrix  of  all  zeroes 

j'k* 

j*,  k*  position. 


except 


1 at  the 


This  approximation  simplifies  the  optimization 
considerably.  However,  this  usually  leads  to  an 
overparameterized  model.  In  other  words,  once  K and  B are 
determined,  it  is  not  possible  to  find  any  corresponding 
p,  Q and  R which  have  the  desired  (known  a priori) 
structure.  A good  estimate  of  elements  in  P # C and  R 
matrices  can  be  obtained  by  a least-squa res  type 
approach.  The  fit  to  the  observed  data  is  better  than  with 
true  values  of  parameters  but  the  parameter  estimates  do  not 
have  minimum  variance.  This  approximation  is  not  good  in 
aircraft  application  vhere  the  structure  of  Q and  R is  known 
fairly  well,  but  is  excellent  where  there  are  many  process 
noise  sources  and  the  characteristics  of  both  Q and  R are 
relatively  unknown  (e.g.,  economic  systems). 

6.  Maximum  Likelihood  With  No  Process  Or  Measurement 
poise 


The  maximum  likelihood  method  can  be  simplified  when 
either  process  noise  or  measurement  noise  are 
absent.  Assuming  no  process  noise  we  can  proceed  as 
follows. 

If  the  process  noise  is  zero  and  initial  states  are 
known  perfectly,  i.e. , w(t)  and  P(0)  are  zero,  the 
covariance  of  the  error  in  the  predicted  state  is  also 
zero.  It  is  clear  from  ( 23  ) that  the  Kalman  gains  are 
zero.  The  innovations  are  the  output  error,  i.e.. 


V"(i)  = y (i)  - h(x(t.)  ,u(t.)  ,e,t.)  ( 60  ) 

ill 
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ao d the  innovation  covariance  is  ( 2t  ) 


B(i)  = R 


( 61  ) 


the  log- likelihood  function  is. 


H 

log  (se(6|z))  = - — V J (i)  R-iV'd)  + log  |R| 

2 1- 1 


( 62  ) 


which  on  optimizing  for  unknown  parameters  in  R yields 


N 

1 r ^ T 

R * Y L V'  (i)  ( 63  ) 

i=1 


The  eguality  in  { 63  ) holds  only  for  those  elements  of  R 
which  are  not  known  a priori.  For  instance,  even  if  B is 
known  to  be  diagonal,  the  right  hand  side  matrix  will  not  be 
diagonal  in  general,  nut  the  off-diagonal  terms  should  be 
ignored  before  they  are  eguated  to  R.  Using  ( 63  ) in  <62  ) 
we  obtain 


N 

log  («C(0|z))  = -~yL  ^ {i)  H_1  y,(i) 

1 i=1 


♦ constant 


< 64  ) 


The  optimizing  function  is  the  same  as  that  for  th  • output 
error  method  except  that  the  measurement  noise  cc variance 
matrix  is  determined  using  ( 63  ) and  is  used  as  the 
weighting  matrix  in  the  criterion  function.  In  the  output 
error  aethod,  the  measurement  noise  is  assumed  known  and  the 
weighting  function  is  arbitrary. 

The  first 
log-likelihood  function 
are 


— “ 3!  (®l  2 

j21og(  (e|z)  ) 

je 

d k 


The  terms  in  the  second  derivative  are  approximated  as 

a*iog<  aeoiz) ) y f av^d)  * t a ✓ <!>  1 

^ aek  ' ‘ ffi  [ a\  R"  aej  J 

( 67  ) 

With  no  measureaent  noise  we  can  proceed  as 
follows.  If  all  states  are  measured  with  no  noise,  the 
covariance  of  the  error  in  state  estimates  is  zero  at  the 


ind  second  derivatives  of  the 
with  respect  to  unknown  parameters 


s>) 


Z/T  . A V 

v w B 

i=  1 


a ■/(i) 


j 


( 65  ) 


N f* 

i-1  L 


H r T 

a v (i>  s_  a v'(i) 

— K"  i 


ae 


aa 


♦ / 


(i)  R- 


1 a*  y/(i)  1 
“ave  kj 


j 


( 66  ) 
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beginning  of  any  time  update. 


P(i-1!i-1)  = 0 


and  x < i—  1 1 i—  1 ) = x(i-l) 


( 68  ) 


It  is  easy  to  show  in  this  case  that  for  fast  sampling  the 
log-likelihood  function  is  quadratic  in  the  difference 
between  measure  values  of  x and  f(x,u,6,t).  The  method 
reduces  to  the  equation  error  method,  the  weight  w being 
chosen  as 


v 


f(x,u,e,t))  (i 


f <x,u,6,t) )Tdt 


( 69  ) 


Thus,  the  maximum  likelihood  method  and  equation-error 
methods  are  equivalent  except  for  the  technique  for  choosing 
the  weighting  matrix. 


F.  AIRCRAFT  PARAMETER  IDENTIFICATION  USING  THE  MAXIMUM 
LIKELIHOOD  METHOD 

In  applying  the  maximum  likelihood  theory  it  is 

necessary  to  consider  the  aircraft  equations  of  mct.j.on  in 

discrete  form.  It  was  assumed  that  the  structure  of  the 

model  was  known.  The  vector  of  unknown  parameters  from  x (0) , 
# 

0*  G,  P , H , Q,  and  R was  denoted  by  0.  The 
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problem,  therefore,  was  to  identify  0 using  the  maximum 
likelihood  method.  In  constructing  the  problem  it  was 
further  assumed  that  the  aircraft  which  had  been  flying  for 
seme  time  was  in  a steady  state  condition,  and  process  noise 
(wind  gust  effects)  was  not  present. 

In  the  computer  implementation,  the  negative 
log-likelihood  function  was  minimized  using  a Gauss-Newton 
gradient  procedure.  A least  sguares  rout' le  was  utilized  to 
provide  approximate  airframe  parameter  start-up 
values.  Additionally  the  results  of  the  SCIDNT  computer 
program  were  utilized  to  obtain  frequency  and  damping. 


G.  SCOPE  OF  RESEARCH 


Research  was  limited  to  the  longitudinal  case  since  both 
available  computer  and  research  time  were  limited.  Although 
additional  data  were  desired,  the  project  airplane  was 
substantially  damaged  by  fire  and  was  removed  from  further 
test  activity.  Further  data  gathering  with  the  F-14A 
airplane  was,  therefore,  not  possible.  Specific 
longitudinal  stability  derivatives  identified  are  given  in 
the  Results  and  Discussion  section.  There  were  no  airplane 
flight  restrictions  which  affected  the  test  methods  or 
limited  the  scope  of  research.  The  research  was  limited  to 
establishing  base  line  data  for  the  short  period  and  phugoid 
modes  and  thence  to  determine  the  feasibility  of  identifying 
these  modes  simultaneously. 


H.  HETHOD  OF  TESTS 

In  general  test  proceedures  were  in  accordance  with 
those  contained  in.  the  0.  S.  Naval  Test  Pilot  School  Fixed 
Wing  Stability  and  Control  Manual,  Ref  5.  There 
were,  however,  some  non-standard  test  techniques  utilized  in 
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exciting  the  desired  airplane  modes.  The  short  period  mode 
was  excited  by  the  standard  doublet  method  and  additionally 
by  applying  a sinusoidal  control  input  at  varying 
frequencies.  The  phugoid  mode  was  excited  by  the  Delta 
Airspeed  method  in  all  cases,  [5]. 

In  exciting  the  phugoid  and  short  period  modes 
simultaneoulsy,  a combination  of  standard  test  techniques 
were  utilized.  The  “simultaneous"  technique  consisted  of 
establishing  the  airplane  in  a phugoid  maneuver  and  during 
the  first  quarter  cycle  of  the  maneuver  applying  a doublet 
input  to  excite  the  short  period  mode. 

The  data  system  utilized  within  the  airplane  was 
specifically  designed  to  be  compatible  with  the  Naval  Air 
Test  Center  Real-Time  Telemetry  Processing  System  (RTPS) . A 
brief  description  of  the  airplane  and  ground  data  systems  is 
presented  in  Appendix  C.  Data  reduction  was  performed  cn  the 
RTPS  Xerox  Sigma  9 computer  in  the  batch  processing 
mode.  Raw  airplane  data  were  converted  to  an  Engineering 
Unit  tape  and  used  as  input  to  the  digital  parameter 
identification  program  (SCIDNT) . In  addition  a least  sguares 
parameter  identification  routine  was  utilized  to  generate 
parameter  start  up  values  for  SCIDNT.  This  least  sguares 
routine  extracted  a set  of  approximate  stability  derivatives 
by  minimizing  the  error  between  the  model  and  actual 
airplane  state  vector  time  histories. 


II.  RESULTS  AND  DISCUSSION 


A.  SHORT  PERIOD  AIRFRAME  PARAMETER  IDENTIFICATION 


1.  Air  plane  Jest  Conditions 


Initial  trimmed  test  conditions  for  each  control 
input  are  presented  in  Table  I.  Prior  to  each  control  input 
the  airplane  was  trimmed  for  steady  state  "hands  off"  flight 


[5J. 


Sine  Nave  Input 
Doublet  Input 


Altitude  MACH 

(ft  HPC) 

29  925  .6149 

30  503  .6289 


CG  Loading 

(%  MAC) 

( i ) 

11.18  B 

10.985  B 


(1) . Pour  Sparrow  III  missiles  (AIM-7)  on  stations  3,  4,  5 

and  6 


Table  I.  Short  Period  Test  conditions 


2.  Method  Of  Investigation 

The  airplane  model  equations  used  were  those  shown 

in  Appendix  B.  In  order  to  isolate  the  short  period  mode 

from  the  longitudinal  equations  of  motion  the  state  vector 

parameters  « , q,  and  8 were  made  insensitive  to  u.  This  was 

accomplished  simply  by  fixing  X^,X,X,Z/H  and  X* 

u q u u e 


46 


equal  to  zeco  and  not  identifying  any  of  then.  This  equation 

modification  resulted  in  the  clasical  short  period 
approximation  Ref.  6. 

To  provide  SCIDNT  with  approximate  airframe 
parameter  start-up  values, the  least  squares  routine  was 
utilized.  For  the  short  period  case  the  least  squares  start 
up  values  were  very  good,  a comparison  of  a typical  set  of 
least  squares  parameters  and  the  final  SCIDNT  parameter 
values  is  shown  in  Table  II. 


Z 

Z 

H 

N 

nr 

oC 

g 

OC. 

q 

Least 

Squares 

-.4909 

0.9624 

-1 . 78 

-.9085 

0.210 

6.09 

SCIDNT 

-.4903 

0.  9805 

-1.95 

-.795 

0.278 

5.542 

Table  II.  Comparison  of  Least  Squares  Parameter 
Values  and  SCIDNT  Parameter  Values 

In  actual  practice,  however,  airplane  data  was 
recorded  in  long  bursts  which  usually  included  three  or  more 
maneuvers  , executed  with  near  identical  initial 
conditions.  As  a result  once  a set  of  start  up  values  were 
determined  for  a given  test  condition  it  was  unnessary  to 
determine  a set  of  start  up  values  for  other  like 
maneuvers. 

The  actual  parameter  identification  proceedure  , for 
the  short  period  mode,  was  relatively  simple  once 
familiarity  with  the  digital  program  was  attained.  Upon 
obtaining  a set  of  "reasonable"  initial  parameter 
estimates,  all  six  short  period  parameters  were  estimated  in 
an  iterative  fashion.  Here  '•reasonable1'  denotes  proper 
parameter  sign  and  magnitude.  Each  computer  run  was  made 
with  the  number  of  iterations  and  step  cuts  both  initialized 
at  ten.  During  a given  iteration  the  program  attempted  to 
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minimize  the  likelihood  function.  If  the  likelihood  function 
increased,  a series  of  step  cuts  would  be  taken  until  the 
likelihood  function  again  decreased  or  the  maximum  number  of 
step  cuts  were  exceeded.  The  program  would  then  continue  to 
iterate  until  the  desired  number  of  iterations/step  cuts 
were  reached.  Due  to  program  logic  one  more  step  cut  than 
reguested  was  executed.  Por  the  short  period 
node,  however,  convergence  was  obtained  early  in  program 
execution.  In  general,  iterations  conducted  with  doublet 
control  input  data  converged  very  quickly  and  began 
executing  step  cuts  until  the  allowable  number  was 
exceeded.  Figure  2 shows  the  change  in  the  likelihood 
function  during  iteration  for  a doublet  input.  Program 
convergence  had  been  essentially  reached  by  the  fourth 
iteration  at  which  time  a series  of  step  cuts  were  made. 


Iterat ion 


Doublet  Input 


The  convergence  characteristics  using  sinusoidal 
input  data  were  somewhat  different.  During  the  program 
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iteration  the  likelihood  function  continuously  decreased  and 
no  step  cuts  were  taken,  Pigure  3.  Program  convergence  for 
this  case  was  evidenced  by  a leveling  out  of  the  likelihood 
function,  but  not  until  the  eighth  iteration.  It  was  felt 
that  since  the  data  length  for  the  sinusoidal  input  was 
longer  than  the  short  period  data  length  the  larger  number 
of  iterations  until  program  convergence  was  not  abnormal. 


Iteration 


Figure  3.  Likelihood  Function  Change  During  Iteration  For 

Sinusoidal  Input 


After  convergence  was  obtained  for  the  parameter 

estimates,  the  Angle  Of  Attack  Scale  Factor  (K^  ) was 

determined.  was  a scaling  factor  in  the  angle  of  attack 

measurement  equation.  Appendix  B.  In  estimating  K two 

terms  in  the  model  were  actually  affected,  itself  and 

the  variable  <1^,  /'J)  which  was  used  to  multiply  g in  the 

measurement  equation.  Since,  in  essence,  two  terms  were 

being  estimated  a modified  form  of  parameter  estimation  was 

utilized.  In  the  term  <1.,  K /0)  the  value  of  K was 

^ ot  x 

initially  set  to  1.0  and  the  single  term  K was 


49 


determined.  This  new  value  of  K was  then  substituted  into 

ct 

the  term  (1^  K /U)  and  again  a set  of  iterations  were 
performed  to  determine  an  updated  . This  aethod  of 

iteration  was  continued  until  there  was  essentially  no 
change  in  , at  which  time  the  value  was  said  to  have 
converged.  In  most  cases  only  two  iterations  were  reguired 
to  obtain  a converged  . 

Upon  determining  a value  for  ,a  final  set  of 

iterations  were  performed  to  determine  values  for 

Z , X , » , 6 . At  this  point,  the  short  period  mode  time 
o o o o 

history  fits  were  nearly  identical  and  the  estimate 
statistics  indicated  a high  degree  of  confidence.  Attempts 
at  further  iterations  to  improve  the  estimates  resulted  in 
using  a lot  more  computer  time , with  very  little 
improvement.  Therefore,  the  procedure  utilized  to  estimate 
the  short  period  mode  was  to  first  estimate  the  short  period 
parameters,  second  estimate  , and  finally,  estimate  the 
initial  values. 

3 • Data  Analy sis 

Data  from  typical  short  period  estimation  runs  for  a 
doublet  and  sinusoidal  input  are  presented  in  Tables  III  and 
IV.  Typical  short  period  time  history  fits  are  presented  in 
Appendix  D Figures  1 and  2,  for  the  doublet  and  sinusoidal 
inputs  respectfully.  It  should  be  noted  that  the  scales  for 
the  short  period  time  histories  are  different  for  the 
doublet  and  sinusoidal  inputs.  Furthermore,  since  the 
estimated  stability  derivatives  were  in  dimensional  form  the 
actual  derivative  values  differed  slightly  from  maneuver  to 
maneuver. 

The  data  obtained  with  the  doublet  input  in  most  all 
cases  converged  to  believable  derivative  estimates,  i.e. , 
expected  sign  and  magnitude,  with  favorable  iteration 
statistics.  For  example,  refering  to  Table  III,  the  value  of 
was  -1.844  with  a standard  deviation  of  0.0061  and  a 
final  step  size  of  -0.0146  which  is  significantly  lower  in 
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magnitude  than  the  estimated  parameter.  Prom  probability 

theory  we  can  say  that  the  probability  is  99.74%  that  the 

value  of  is  -1.844  ±0.0183  vhich  is  veil  within 

cC 

acceptable  engineering  tolerances.  This  same  analysis  was 
carried  out  for  the  other  estimated  derivatives  with  similar 
results. 


Derivatives 

Z 

eC 

Z 

q 

M 

oC 

n 

g 

Estimate 

-.4138 

1.099 

-1.844 

-.826 

0.1628 

5.098 

Standard 

Deviation 

0.0024 

0. 0038 

0. 0061 

0.0037 

0.0055 

0.0104 

Parameter 
Step  Size 

-.0006 

-.012 

-.0146 

-.0001 

0.0143 

0.028 

Eigenvalue 

Spread 

1.28 

x 10* 

Likelihood 

Function 

Value 

-5.  1 

X 103 

Time 

History 

Matches 

o 

N 

6* 

.00396 

it 

0.564 

<r  * o. 

9 

0039 

ii 

0.0043 

0~  ~ 2. 
n 

z 

33 

Table  III.  SCIDNT  Estimates  of  Longitudinal  Short  Period 
Stability  Derivatives  Using  a Doublet  Input 


The  eigenvalue  spread,  that  is,  the  difference 
between  the  largest  eigenvalue  and  the  smallest  eigenvalue 
of  the  information  matrix,  was  1.28  x 10*  vhich  was  less 
than  the  maximum  spread  allowed  by  SCIDNT  (10*°).  All 
parameters  converged  indicating  no  identif iabilit y 
problems,  and  the  eigenvalue  spread  was  therefore  not 
considered  significant.  The  likelihood  function  value  was  a 
large  negative  number  which  converged  rapidly  (Pigure 
2),  further  indicating  no  identif iability  problems. 
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As  a final  check  of  stability  derivative  correctness 
a series  of  tine  history  fits  were  constructed.  In 
constructing  these  plots  both  the  model  and  airplane  were 
excited  with  the  same  control  input  and  the  airplane  time 
history  plots  superimposed  on  the  same  scale.  The  time 
history  fit  for  the  doublet  input  was,  in  general,  very 
good,  with  the  difference  between  the  measured  value  and 
estimated  value  within  the  instrumentation  measurement  error 
tolerances  Ref.  7. 

Referring  to  Table  IV,  it  can  be  seen  that  typically 
the  estimated  derivatives  using  a sinusoidal  input  converged 
to  appropriate  (sign  and  magnitude)  stability  derivative 
values. 


Derivatives 

Z 

<£ 

z 

g 

M 

cc 

t! 

g ' e 

Estimate 

-.4303 

1. 1229 

-2*007 

-.7615  0.1532 

5.2274 

Standard 

Deviation 

0.0042 

0.0061 

0.013 

0.007  0.0045 

0.0131 

Parameter 
Step  Size 

0.0007 

-.0016 

-.0043 

-.0061  -.0037 

0.0332 

Eigenvalue 

Spread 

1.67 

x IQ* 

Likelihood 

Function 

Value 

-7.323  x103 

Time 

History 

Matches 

O'  - 0 

tX 

.00533 

cr  = 
u 

0.9  <T  = 0 . 

q 

006 

n 

0.01 

T = 2.32 
n 

z 


Table  IV.  SCIDNT  Estimates  of  Longitudinal  short  Period 
Stability  Deravitives  Using  A Sinusoidal  Input 


In  general,  however,  the  parameter  standard  deviations  and 
step  sizes  were  not  as  favorable  as  those  obtained  with  the 
doublet  input  for  these  two  examples.  It  was 

52 


possible,  however,  to  select  another  pair  of  maneuvers  and 
reach  the  opposite  conclusion,  ie.  sinusoidal  input  more 
favorable.  Since  all  parameter  estimates  for  both  the 
doublet  and  sinusoidal  inputs  converged  it  was  concluded 
that  the  slight  difference  in  statistical  data  was  not 
significant.  If  a determination  as  to  which  control  input  to 
use  were  required,  it  was  felt  that  the  primary 
consideration  would  be  ease  of  pilot  input.  If  an  autopilot 
device  were  available  which  could  input  the  sinusoid  then 
this  input  should  be  used,  however,  with  no  autopilot  system 
installed  the  doublet  input  would  provide  excellent 
results. 

A consistent  problem  that  occured  throughout  the 

research  was  concerned  with  the  initial  conditions  prior  to 

the  control  input.  If  the  initial  condition  of  the  airplane 

was  not  stabilized  prior  to  the  control  input  or  if  a data 

burst  was  taken  in  the  middle  of  a control  input,  no  time 

history  match  was  possible.  According  to  theory  it  is 

possible  to  extract  the  airframe  parameters  at  any  time 

during  a maneuver,  this,  however  did  not  appear  to  be 

true.  Fi.  ire  3 Appendix  D shows  the  results  of  estimating 

the  parameters  starting  from  non-zero  initial 

conditions.  This  time  history  fit  was  generated  after  all 

iterations  on  both  parameters  and  initial  conditions  were 

completed.  The  parameter  estimates  did  not  converge  and  the 

most  powerful  derivative,  M , was  much  larger  in  magnitude 

c t- 

than  any  ether  estimates  with  a value  of  -3.7.  These  results 
occured  any  time  an  attempt  was  made  to  start  SCIDNT  from  an 
initial  condition  other  than  steady  state.  It  was  concluded 
that,  within  the  scope  of  this  research,  meaningful 
parameter  estimates  were  only  obtainable  when  started  from 
an  initial  steady  state  condition. 

A comparison  of  short  period  characteristics 
obtained  from  Navy  Preliminary  Evaluation  (NPE)  data  and 
those  determined  by  SCIDNT  is  presented  in  Table  V. 
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Description 

Damping 

Ratio 

Frequency 

t 

07 

bsp 

n 

sp 

(rad/sec) 

NPE 

0.35 

1.49 

SCIENT-Doublet 

0.403 

1.54 

SCIDNT-Sinusoidal 

0.371 

1.61 

Table  V.  Comparison  of  Longitudinal 
Short  Period  Characteristics 


It  was  concluded  that,  within  the  scope  of  this 
research,  the  parameter  identification  program  SCIDNT  was 
easy  and  straightforward  to  use  in  determining  the  short 
period  mode.  In  addition,  it  was  concluded  that  it  was  not 
possible  to  obtain  accurate  estimates  when  starting  from 
non-steady  state  initial  conditions. 


B.  PHUG01D  AIRFRAME  PARAMETER  IDENTIFICATION 

1 • Airplane  Test  Conditions 

Initial  trimmed  test  conditions  for  the  phugoid 
maneuver  are  presented  in  Table  VI.  Prior  to  the  control 
input  the  airplane  was  trimmed  for  steady  state  "hands  off" 
flight  [ 5 ]. 


Altitude  MACH 

CG 

Loading 

(Ft  HPC ) 

(%  MAC) 

Phugoid  Input  30  055  0.614 

10.94 

B 

Table  VI.  Phugoid  Test 

Conditions 
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2.  Method  Of  Investigation 


The  airplane  model  equations  used  were  those  shown 

in  Appendix  B.  The  derivitives  X and  X^,  were  fixed  at 

oc  ve 

zero  values  since  their  effect  on  the  phugoid  was 

small.  Additionally,  it  was  felt,  fixing  these  terms  would 

facilitate  the  identification  of  the  dominant  phugoid 

terms.  Finally  from  Appendix  B it  can  be  seen  that  the  tern 

X'  was  defined  as  (X  * oc  l)  ) and  therefore  due  to  its 
q q O o 

relatively  large  size  was  approximated  by  cC  0 and 

o o 

fixed.  The  parameters  which  were  identified  during  the  pure 

phugoid  maneuver  were  therfore,  Z , X , and  M . 

u u u 

Since  there  was  no  short  period  input  during  the 
pure  phugoid  maneuver  it  was  not  possible  to  effectively 
identify  the  short  period  parameters.  Additionally  the 
program  solution  diverged  without  a set  of  predefined  short 
period  parameters  and  it  was  necessary  to  obtain  parameters 
from  a prior  run.  A set  of  short  period  parameters 

were,  therefore,  obtained  from  a previous  identification  run 
and  fixed  throughout  the  phugoid  parameter  identification 
runs. 

In  an  attempt  to  provide  SCIDNT  with  approximate 

airframe  parameter  start-up  values  the  least  squares  routine 

was  utilized.  For  the  phugoid  case  the  least  sguares  start 

up  values  were  unusable.  For  example,  least  squares  values 

for  X_  was  -2.6697  and  XA  was  -240.5,  both  of  which  are 
* °e 

terms  of  small  absolute  value.  As  an  additional  check,  the 

least  squares  parameter  estimates  were  used  to  generate  a 
set  of  time  histories.  For  these  time  bistorys  the  standard 
deviation  of  airspeed  and  normal  acceleration  was  49.12 
ft/sec  and  5.44  ft/sec2  respectively.  Therefore,  since  the 
least  sguares  routine  did  not  provide  acceptable  start-up 
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values  the  identified  parameters  were  set  initally  to  zero 
for  the  first  SCIDNT  iteration.  Further  attempts  to  obtain 
phugoid  start  up  values  for  SCIDNT  from  the  least  squares 
routine  proved  futile  and  its  use  was  discontinued  for  the 
phugoid  case. 

The  version  of  SCIDNT  utilized  for  this 

investigation  would  accept  a maximum  of  500  data  samples  per 
run.  Since  the  phugoid  period  was  approximately  100  sec  a * 
data  burst  recorded  at  20  samples  per  second  (SPS)  would 
only  include  the  first  25  sec  of  the  maneuver.  This  time 
span  was  obviously  acceptable  for  the  short  period  mode  but 
clearly  unacceptable  during  a phugoid  maneuver.  An  attempt 
was  therefore,  made  to  utilize  data  at  2 and  5 SPS  to 
include  as  much  of  the  phugoid  maneuver  as  possible.  This 
data  was  from  the  same  maneuver  with  data  tapes  recorded  at 
2 and  5 SPS. 

Although  the  2 SPS  data  would  provide  a longer 
phugoid  time  spar,  in  theory,  the  program  exhibited  severe 
instabilities  during  the  control  input  and  therefore  only 
the  latter  portion  of  the  data  burst  was  usable.  It  was 
strongly  suspected  that  these  instabilities  were  due  to  the 
unusually  rough  control  inputs. 

The  actual  parameter  identification  procedure  for 
the  phugoid  mode  was  somewhat  different  than  that  utilized 
during  the  short  period  investigation.  The  short  period 
stability  derivatives  were  obtained  from  a run  whose  initial 
conditions  most  closely  approximated  those  of  the  phugoid 
run.  These  short  period  parameters  were  then  fixed  as 
constants  during  the  phugoid  iterations.  Computer  runs  were 
made  with  the  number  of  iterations  and  step  cuts  both 
initialized  at  ten.  The  program  would  then  continue  to 
iterate  until  the  desired  number  of  iterations/step  cuts 
were  completed.  For  both  the  5 SPS  and  2 SPS  data 
convergence  was  obtained  prior  to  executing  the  full  10 
iterations.  Iterations  conducted  with  phugoid  data  recorded 
at  2 SPS  appeared  to  converge  very  guickly  and  began 
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executing  step  cuts  until  the  allowable  nunber  was 
exceeded.  As  in  the  short  period  case  the  program  executed 
one  step  cut  more  than  requested.  Pigure  4.  presents  the 
change  in  the  likelihood  function  during  iteration  using 
phugoid  data  recorded  at  2 SPS.  The  program  converged  to  its 
"best"  solution  by  the  third  iteration  at  which  time  the 
step  cuts  were  executed. 


Iteration 


Pigure  4.  Likelihood  Function  During  Iteration 

Using  2 SPS 


Note  that  although  the  program  had  apparently 
converged  at  the  third  iteration  the  value  of  the  likelihood 
function  decreased  until  the  fourth  step  cut.  The  function 
value  then  remained  essentially  constant  throughout  the 
remaining  step  cuts. 

The  convergence  characteristics  using  data  recorded 
at  5 SPS  were  somewhat  different.  Figure  5. 
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Iteration 


Figure  5.  Likelihood  Function  During  Iteration 

Using  5 sps 


It  should  be  noted,  however.,  that  the  scales  are 

different  for  each  plot.  The  program  had  essentially 

converged  at  the  fifth  iteration,  at  which  time  the  program 

executed  step  cuts.  Note  that  the  likelihood  function 

values  decrease  using  5 SPS  data,  and  were  essentially  the 

same  using  the  2 SPS  data.  However,  the  time  history  matches 

for  each  set  of  data  were  significantly  different.  The  angle 

of  attack  scale  factor  (K  ) was  not  estimated  since  the 

cc 

value  of  from  the  short  period  data  was  utilized.  After 

convergence  was  obtained  for  the  parameter  estimates  a final 
set  of  iterations  was  performed  to  estimate  values  for 

Z , X , H , 0 and  the  three  phi'goid  parameters 


simultaneously . At  this  point,  the  phugoid  tine  history  fits 

mere  acceptable  for  the  5 SPS  data  and  clearly  unacceptable 
using  the  2 SPS  data.  Further  attenpts  to  inprove  the 
estimates  and  tine  histories  resulted  in  using  a lot  more 
computer  tine  with  no  improvement. 

3.  pat?  Analysis 

Data  fron  phugoid  estimation  runs  using  2 SPS  and  5 

SPS  are  presented  in  Tables  VII  and  VIII 

respectfully.  Phugoid  time  history  fits  are  presented  in 

Appendix  D Figures  4 and  5 for  2 and  5 SPS  respect  fully,.  It 

should  be  noted  that  the  scales  of  the  phugoid  time 

histories  for  2 and  5 SPS  are  different. 

The  data  obtained  using  a sample  rate  of  2 SPS  was 

totally  unacceptable.  For  example,  refering  to  Table 

VII,  the  value  of  M was  -0.00000017  while  the  last  step 

u 

size  was  0.003. 


Derivatives 

Z 

u 

X 

cC 

X 

u 

X 

q 

u °e 

Estimate 

-.003 

0 

-.758 

-47.9 

-1.7x10-*  0 

Standard 

deviation 

.00006 

fixed 

.0542 

fixed 

.0002  fixed 

Parameter 
Step  Size 

-.001 

fixed 

0.444 

fixed 

0.0C3  fixed 

Eigenvalue 

Spread 

3.54 

x 10* 

Likelihood 

Function 

Value 

— 8.4 

X 103 

Time 

History 

Matches 

ii 

o 

.004  8 

<r  = 

u 

26.77 

<r  = 0.00534 

q 

n 

0.053 

. 

(N 

II 

53 

z 


Table  VII.  SCIDNT  Estimates  of  Longitudinal  Phugoid 
Stability  Derivatives  Using  2 SPS  Data 


Although  the  eigenvalue  spread  and  likelihood 
function  change  during  iteration  indicate  no  significant 
program  difficulties  the  standard  deviations  of  the  state 
error  vector  are  very  large.  As  an  approximation  it  can  be 
said  that  the  phugoid  maneuver  is  essentially  a varying 
airspeed  and  pitch  attitude  maneuver  with  a constant  load 
factor.  From  Table  VII  it  can  be  seen  that  the  standard 
deviation  for  airspeed  was  26.77  ft/sec-1  and  for  pitch 
attitude  was  0.053  rad,  both  of  which  were  unacceptable. 

The  time  history  plots  constructed  using  the  2 SPS 
data  showed  no  real  agreement  between  the  model  and  the 
airplane.  In  constructing  these  plots  both  the  model  and 
airplane  were  excited  with  the  same  control  input  and  the 
airplane  time  history  plots  superimposed  on  the  same 
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scale.  It  is  obvious  that  there  was  simply  no  fit  of  any  of 
the  state  vector.  It  should  be  noted  that  the  plot  of 
average  horizontal  stabilizer  position  although  appearing  to 
be  very  rough  was  a result  of  the  plot  scale  only  and  was  in 
essense  relatively  smooth. 

The  parameter  identification  results  using  the  5 SPS 
data  indicated  it  was  possible  to  identify  the  phugoid  mode 


parameters. 

Referring 

to  Table  VIII 

, it  can  be  seen 

that  the 

estiaat  ed 

derivatives  using  the 

5 SPS  data  converged  to 

appropriate 

(sign  and  magnitude} 

stability  derivative 

values. 

Der ivatives 

Z 

u 

X 

R 

X 

c 

x a 

q U 

Estimate 

-.0002 

0 -.006 

-45.5  -9.6x10 

-5  0 

Standard 

Deviation 

1.2x10-* 

fixed  .0002 

fixed  5.x10-6 

fixed 

Parameter 
Step  Size 

-8.7x10-6 

fixed  -.0016 

fixed  -.00002 

fixed 

Eigenvalue 

Spread 

1.82  x 10*2 

Likelihood 

Function 

Value 

-1. 

49  x 10* 

Time 

History 

Hatches 

• 

o 

II 

b* 

00267  (T 

u 

= 5.53  Q-  = 0. 

9 

00343 

<T  = 0.014 
6 

CT  = 0.92 
n 

z 


Table  VIII.  SCIDNT  Estimates  of  Longitudinal  Phugoid 
Stability  Derivatives  using  5 SPS  Data 


The  standard  deviations  of  the  parameter  estimates  were  much 
smaller  than  the  parameter,  and  the  parameter  step  sizes 
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indicated  that  the  estimates  had  converged  to  the  "best’' 
estimates  obtainable  with  the  data.  It  should  be  noted  that 
standard  deviations  for  u and  8 were  5.53  ft  sec-1  and  0.014 
rad  respectively/  a significant  reduction  in  standard 
deviation  from  those  obtained  with  2 S PS  data 

The  tine  history  plots  using  5 SPS  data  showed 
relatively  close  agreement  between  the  airplane  and  model 
time  histories.  There  was,  however,  some  model  instability 
noted  during  the  first  30  seconds  of  the  time  histories.  It 
was  concluded  that  these  instabilities  were  due  to  the 
relatively  slow  sampling  rate  during  the  initial  control 
input. 

From  SCIDNT  the  phugoid  damping  ratio  ( f ) and 

bp 

frequency  (co  (rad  sec-1))  were  0.0486  and  0.06 
n 

respectfully.  Available  P-14  data  was  obtained  in  the  auto 

sweep  mode  with  SAS  on  while  SCIDNT  data  was  obtained  in 
auto  sweep,  SAS  off.  It  was  therefore  not  possible  to  make 
any  meaningful  comparisons.  However,  it  was  noted  that 
between  10,000  ft  and  35,000  ft  the  phugoid  period  varied 
between  106  sec  and  72  sec.  The  SCIDNT  frequency  resulted  in 
a period  of  92.53  sec  which  when  considered  in  light  of 
available  N PE  data  appeared  reasonable. 

It  was  concluded  that,  within  the  scope  of  this 
research,  the  paramater  identification  program  SCIDNT 
provided  a reasonable  estimate  of  the  phugoid  longitudinal 
parameters. 


C.  SIMULTANEOUS  SHORT  PERIOD  AND  PHUGOID  PARAMETER 
IDENTIFICATION 

1.  Airplane  Test  Conditions 

Initial  trimmed  test  conditions  for  the  short  period 
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and  phugoid  maneuvers  are  presented  in  Table  IX.  Prior  to 
the  control  input  the  airplane  was  trimmed  for  steady  state 
"hands  off"  flight  [5],  Since  ths  short  period  input  was 
applied  after  the  phugoid  input,  the  two  test  conditions  are 
different. 


Altitude 
(Ft  HPC ) 

HACH 

CG 

(X  MAC) 

Loading 

Phugoid  Input 

30  326 

.6217 

10.86 

B 

Short  Period  Input 

30  703 

.5963 

10.86 

B 

Table  IX.  Short  Period  and  Phugoid  Test  Conditions 

2.  Method  Of  Investigation 

The  airplane  model  equations  used  were  those  shown 

in  Appendix  B.  The  derivatives  X , , and  X were  fixed 

oc  °e  g 

as  described  in  the  previous  section. 

The  maneuver  investigated  in  this  section  was 
designed  specifically  to  provide  the  means  of  identifying 
the  short  period  and  phugoid  parameters  from  one  maneuver 
or,  if  you  will,  simultaneously.  Initiation  of  the  phugoid 
maneuver  was  accomplished  in  accordance  with  standard  test 
techniques,  [5].  Upon  establishing  the  airplane  in  the 
phugoid  maneuver,  a doublet  input  was  applied  to  excite  the 
short  period  mode.  The  airplane  was  then  allowed  to  fly 
through  at  least  one  phugoid  cycle. 

Identification  of  the  short  period  parameters  was 
essentially  routine  and  straightforward.  The  procseaures 
used  in  estimating  the  parameters  were  those  described  in 
the  previous  short  period  section.  The  computer  runs  were 
made  with  the  number  of  iterations  and  step  cuts  both 
initialized  at  ten.  Although  this  short  period  input  was 
superimposed  on  the  phugoid,  convergence  was  obtained  early 
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in  program  execution.  Iterations  conducted  under  these 
conditions  converged  very  quickly  and  executed  step  cuts 
until  the  allowable  number  was  exceeded.  Likelihood  function 
change  during  the  short  period  iterations  was  similar  to 
that  shown  in  Pigure  2. 

Based  on  the  results  of  the  previous  section  it  was 
concluded  that  data  recorded  at  5 SPS  would  most  likely 
yield  reasonable  phugoid  parameter  estimates.  The  short 
period  stability  derivatives  obtained  from  this  run  were 
fixed  as  constants  during  the  phugoid  iterations.  Computer 
runs  were  then  made  with  the  number  of  iterations  and  step 
cuts  both  initialized  at  ten. 

Figure  6 presents  the  likelihood  function  value 
during  iteration  of  the  phugoid  stability  derivatives. 


I teration 


Pigure  6.  Likelihood  Function  During  Iteration 

For  Phugoid 
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The  likelihood  function  change  for  the  phugoid  maneuver 
displayed  some  characteristics  heretofore  not 

encountered.  The  characteristics  shown  are  typical  of  the 
results  obtained  with  the  available  data,  At  some  time 
during  the  iteration  sequence  the  likelihood  function  would 
exhibit  a jump  discontinuity,  in  this  case  it  occurred 
between  the  first  and  second  iterations.  Program  iteration 
then  appeared  to  proceed  in  a normal  manner  until  the  tenth 
iteration  at  which  time  the  program  terminated.  From  Figure 
6 it  was  obvious  that  convergence  had  not  been  obtained  at 
program  termination. 

In  an  attempt  to  correct  the  iteration  instabilities 

the  phugoid  stability  derivatives  2 , X , and  Fl  were 

u u u 

estimated  simultaneously  with  2 , X , ti  , and  0 . Figure  7 

o o o o ' 

presents  the  likelihood  function  during  a typical  iteration 
for  seven  parameters. 
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Iteration 


Figure  7,  likelihood  Function  During  Iteration  lor 

Pbugoia  (Stability  Derivatives,  Z , X , rt  , and  e ) 

o o o o 

It  can  be  seen  that  the  likelihood  function  value  converged 
rapidly  to  a steady  state  value  with  no  apparent 
discontinuities  or  nonlinearaties. 

3.  Data  Analysis 

Data  from  a typical  short  period  run  superimposed  on 
the  phugoid  are  presented  in  Table  X.  A typical  time  history 
fit  is  presented  in  Appendix  D,  Figure  6. 
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Derivatives 

Z 

at 

Z 

q 

M 

on 

H 

g 

*e 

% 

Estimate 

-.4285 

1.048 

-1.54 

» . 665 

0. 19 

4.464 

Standard 

Deviation 

0.0049 

0.0074 

0.014 

0.011 

0.0082 

0.0312 

Parameter 
Step  Size 

0.0149 

0.0119 

-.0166 

0.0355 

0.0173 

-.1815 

Eigenvalue 

Spread 

1.25 

X 105 

Likelihood 

Function 

Value 

-4.36  x 103 

Time 

History 

Matches 

o 

ii 

b« 

.00188 

II 

0.8822 

. 

o 

it 

b01 

0036 

O'  = 

6 

0.0048 

O'  * 1. 
n 

82 

z 


Table  X.  SCIDNT  Estimates  of  Longitudinal  Short  Period 

Stability  Derivatives 


The  short  period  parameter  estimates  displayed 
characteristics  similar  to  those  observed  when  estimating 
the  short  period  parameters  from  a short  period  input 
only.  The  time  history  plots  shoved  good  agreement  between 
the  model  and  actual  airplane  data.  It  was  noted  that 
airspeed  was  not  constant  since  the  short  period  was 
superimposed  on  the  phugoid.  This  varying  airspeed  did  not 
however,  have  any  noticable  effect  on  the  determination  of 
the  short  period  parameters. 

The  results  for  the  phugoid  portion  of  the  maneuver 
were  at  times  in  disagreement  with  theory.  From  theory  it 
can  be  shown  that  the  speed  damping  derivative,  which  acts 
along  the  negative  X axis  increases  with  speed  and  is 
therefore  negative  [6],  During  the  iteration  to  determine 
the  phugoid  stability  derivatives  it  was  noted  that  the 


67 


' . I 

i 

» 

’ speed  damping  derivative  (X  ) , became  positive  on  the 

seventh  iteration  and  remained  positive  until  program 


termination.  Table  XI. 


Derivatives  Z 

u 

X 

Ct 

X 

u 

x a 

q u ve 

Estimate  -.0002 

0 

.005 

-43.85  -7.8x10-*  0 

Standard  1.1x10“* 

Deviation 

fixed 

.0002 

fixed  3.5x10“*  fixed 

Parameter  5.38x10“* 
Step  Size 

fixed 

4.8x10“ 

s fixed  1.2x10-s  fixed 

Eigenvalue 

Spread 

6.2 

x 10>2 

Likelihood 

Function 

Value 

-1. 

005  x 10* 

Time 

History  C - 0 

Hatches  cc 

.004 

c r 

u 

= 1.97  O'  = 0.004  7 

q 

<r 

= 0.031 

cr  = i.48 

9 n 

z 

Table  XI.  SCIDNT  Estimates  of  Phugoid  Stability 

Derivatives  (Positive  X ) 

u 

The  time  history  plots  were  in  reasonable 

agreement,  Appendix  D Figure  7.  The  positive  value  of  was 

not  however  considered  reasonable  since  this  value  would 

imply  decreasing  drag  with  increasing  airspeed.  The  X term 

u 

was  however  negative  through  the  sixth  program  iteration.  As 

mentioned  previously  the  likelihood  function  during 
iteration  showed  definite  signs  of  instability,  it  was 
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therefore  suspected  that  the  speed  damping  sign  change  was 
due  to  this  program  instability.  The  program  was 
therefore,  run  again  and  terminated  at  the  sixth  iteration 
thereby  producing  a negative  speed  damping  term.  Table  XII. 


Derivatives  Z 

u 

X 

cc 

X X 

u g 

a 

u 

Estimate  -.00008 

0 

-.00032  -43.85 

0.00017 

0 

Standard  2.6x10-* 

Deviation 

fixed 

0.00047  fixed 

6.7x10-* 

fixed 

Parameter  4.0x10-5 
Step  Size 

fixed 

-.0059  fixed 

0.000058 

fixed 

Eigenvalue 

Spread 

8.53  x 10»i 

Livelihood 

Function 

Value 

-7.68  x 103 

Time 

History  0~  - 0. 

Matches  OC 

0071 

0~  - 3.94 
u 

cr  =0.0053 

g 

o~ 

e 

= 0.0366  0~  -2. 
n 

14 

z 

Table  XII.  SCIDNT  Estimates  of  Phugoid  stability 

Derivatives  (Negative  X ) 

u 


With  the  iterations  thus  restricted  to  six  the  value  for  X 

u 

was  maintained  at  a negative  value.  Time  history  Flots  for 

the  negative  X case  are  presented  in  Appendix  D,  Figure 
u 

8.  It  was  noted  that  there  was  no  appreciable  improvement  in 

the  time  history  plots.  In  fact  from  Tables  XI  and  XII  it 
can  be  seen  that  the  standard  deviation  of  the  state  error 
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case.  This 


vectors  was  degraded  for  the  negative  X 

u 

condition  in  fact  is  however#  not  very  startling  when  it  is 

realized  that  the  program  had  demonstrated  instabilities  and 
had  not  converged. 

In  an  effort  to  reconcile  these  inconsistencies  and 

lack  of  convergence  the  three  phugoid  parameters  were 

estimated  along  with  2 # X , M # and  6 . This  technigue 

o o o o 

yielded  reasonable  stability  derivatives  and  acceptable  time 

history  plots.  Stability  derivative  data  provided  by  the 
seven  parameter  identification  are  shown  in  Table  XIII#  time 
history  plots  are  presented  in  Appendix  D,  Figure  S. 


Derivatives  Z X 

u oc 

X 

u 

X 

q 

M X^ 

u ^e 

Estimate  -.00019  0 

-.0029 

-43.85  -. 

00014  1 0 

Standard  1.4x10-*  fixed 

Deviation 

.00026 

fixed  5. 

2x10“*  fixed 

Parameter  1.23x10~7  fixed 
Step  Size 

6.3x10- 

s fixed  -4. 

02x10~6  fixed 

Eigenvalue 

Spread 

3.7 

x 10*  3 

Likelihood 

Punction 

Value 

-4. 

23  X 103 

Time 

History  O'  = 0.0045 

Matches  cc 

O' 

u 

= 4.17  cr 

<3 

=0.00523 

O' 

e 

= 0.031 

O'  = 1.86 
n 

z 

Table  XIII.  SCIDNT 

Estimates 

of  Phugoid 

Stabili ty 

Derivatives  (Parameters  Estimated  Bith  Z , x # M # and  6 ) 

o o o o 
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The  parameter  standard  deviations  were  small  and 
within  engineering  tolerance.  The  parameter  step  sizes  were 
gaite  small  and  shoved  further  evidence  of  solution 
convergence.  The  standard  deviation  of  the  calculated  model 
response  with  respect  to  the  measured  airplane  response  was 
acceptable.  The  time  history  plots  provided  further 
indication  of  parameter  estimate  validity. 

The  problems  associated  with  the  phugoid  using  the 
simultaneous  technique  could  have  been  a result  of  program 
deficiencies  or  airplane  anomalies  or  a combination  of 
both.  The  flight  conditions  at  which  the  data  was  obtained 
were  areas  in  which  the  phugoid  mode  was  lightly  damped  or 
neutrally  damped.  In  addition  the  airplane  wing  was  set  to 
the  auto  sweep  mode.  It  was  felt  that  either  of  these 
conditions  could  have  caused  the  program  instabilities 
noted.  Attempts  to  obtain  new  data  were  not  possible  due  to 
airplane  non-availability  as  mentioned  earlier. 

It  was  concluded  that  extraction  of  the  short  period 

stability  parameters  in  the  presence  of  the  phugoid  was 

easy,  straight  forward,  and  yielded  results  similar  to  those 

obtained  from  pure  short  period  data.  Estimation  of  the 

phugoid  stability  parameters  was  possible  when  they  were 

estimated  in  conjunction  with  Z , X , H , and  0 . 

o o o o 

It  was  recommended  that  a new  set  of  data  be 
obtained  at  a flight  condition  where  the  phugoid  is  at  least 
moderately  damped  and  that  this  data  be  analyzed  to  resolve 
the  present  anomalies. 


71 


III.  CONC’ USIONS 


A.  GENERAL 

The  maximum  likelihood  identification  program  SCIDNT 
showed  great  promise  in  obtaining  the  short  period  and 
phugoid  stability  derivatives  from  one  maneuver, 
••simultaneously". 

B.  SPECIFIC 

Within  the  scope  of  this  research,  the  parameter 
identification  program  SCIDNT  was  easy  and  straight-f crward 
to  use  in  determining  the  short  period  mode. 

Meaningful  parameter  estimates  were  obtainable  only  when 
started  from  an  initial  steady  state  condition. 

Within  the  scope  of  this  research,  the  parameter 
identification  program  SCIDUT  provided  a reasonable  estimate 
of  the  phugoid  longitudinal  parameters. 

Extraction  of  the  short  period  stability  parameters  in 
t.he  presence  of  the  phugoid  was  easy,  straight-forward,  and 
yielded  results  similar  to  those  obtained  from  pure  short 
period  data. 

Estimation  of  the  phugoid  stability  parameters  was 
possible  when  they  were  estimated  in  conjunction  with 

2 , X , M , and  6 . 
o o o o 
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IV.  RECOMMENDATIONS 


It  is  recommended  that  a new  set  of  data  be  obtained  at 
a flight  condition  where  the  phugoid  is  at  least  moderately 
damped  and  that  these  data  be  analyzed  to  resolve  the 
present  anomalies. 
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APPENDIX  A 

TEST  AIRCRAFT  DESCRIPTION 


The  F-14  airplane  is  a supersonic,  two-place, 
twin-engine,  swing-wing,  air-superiority  fighter  designed 
and  manufactured  by  Grumman  Aerospace  Corporation.  In 
addition  to  its  primary  fighter  role,  carrying  missiles  and 
an  internal  20-millimeter  gun,  the  airplane  is  designed  for 
fleet  air  defense  and  ground  attack  missions.  Photographs  of 
the  test  airplane  are  shown  in  Figures  1 and  2.  A three-view 
drawing  is  presented  in  Figure  3. 

The  forward  fuselage,  containing  the  crew  and  electronic 
equipment,  projects  forward  frOB  the  mid  fuselage  and  wing 
glove.  Outboard  pivots  in  the  highly  swept  wing  glove 
support  the  movable  wing  panels,  which  incorporate  integral 
fuel  cells  and  full-span  leading  edge  slats  and  trailing 
edge  flaps  for  supplemental  lift  control.  In  flight,  the 
wings  may  be  varied  in  sweep,  area,  camber,  and  aspect  ratio 
by  selection  of  any  leading-edge  sweep  angle  between  20°  and 
68».  Hing  sweep  can  be  automatically  or  manually  controlled 
to  optimize  performance  and  thereby  enhance  aircraft 
versatility.  Separate  variable-geometry  air  inlets,  offset 
from  the  fuselage  in  the  glove,  direct  primary  airflow  to 
two  TF30-P-412A,  dual-axial-com pressor , turbo-fan  engines 
equipped  with  afterburners  for  thrust  augmentation.  The 
displaced  engine  nacelles  extend  rearward  to  the  tail 
section,  supporting  the  twin  vertical  tails,  horizontal 
tails,  and  ventral  fins.  The  middle  and  aft  fuselage,  which 
contains  the  main  fuel  cells,  tapers  off  in  depth  to  the 
rear  where  it  accommodates  the  speed  brake  surface  and 
arresting  hook.  Retractable  vanes  in  the  glove  leading  edges 


extend  to  supplement  lift  and  compensate  for  changes  in  the 
aircraft  aerodynamic  center.  All  control  surfaces  are 
positioned  by  irreversible  hydraulic  actuators  to  provide 
desired  control  effectiveness  throughout  the  flight 
envelope.  Stability  augmentation  features  in  the  flight 
control  system  enhance  flight  characteristics  and  provide  a 
more  stable  and  maneuverable  weapons  delivery  platform.  The 
tricycle-type,  forward-retracting  landing  gear  is  designed 
for  nose  gear  catapult  launch  and  carrier  landings.  Missiles 
and  external  stores  are  carried  from  eight  hardpoint 
stations  on  the  center  fuselage  between  the  nacelles  and 
under  the  nacelles  and  wing  glove;  no  stores  are  carried  on 
the  movable  portion  of  the  wing.  The  fuel  system 
incorporates  both  inflight  and  single-point  ground  refueling 
capabilities. 

A more  complete  description  of  the  F-14A  airplane  is 
given  in  Ref.  8. 
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APPENDIX  B 


DERIVATION  OP  THE  AIRPLANE  MODEL  AND  MEASUREMENT  EQUATIONS 

A detailed  derivation  of  the  complete  aircraft  equations 
of  notion  is  shown  in  Ref.  6.  A summary  of  the  assumptions 
made  during  that  derivation  is  shown  below  for 
completeness.  These  assumptions  lead  to  a set  of  force  and 
moment  eguations. 

* Summary  of  Derivation  Assumptions  * 

1.  The  airframe  is  a rigid  body 

2.  The  earth  is  fixed  in  space 

3.  The  earth's  atmosphere  is  fixed  with  respect  to  the 
earth 

4.  The  mass  of  the  airplane  remained  constant  for  the 
duration  of  any  particular  dynamic  analysis 

5.  The  xz  plane  is  a plane  of  symmetry 

6.  The  disturbances  from  the  steady  flight  condition  is 
small  enough  so  that  the  products  and  squares  of  the  changes 
in  velocities  are  negligible  in  comparison  with  the  changes 
themselves 

7.  The  disturbance  angles  are  small  enough  such  that  the 
sines  of  these  angles  may  be  set  equal  to  the  angles  and  the 
cosines  set  equal  to  one.  Products  of  these  angles  are  also 
approximately  zero  and  can  be  neglected 

8.  Since  the  disturbances  were  small,  the  change  in  air 
density  encountered  by  the  airplane  during  any  disturbance 
can  be  considered  zero 
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9.  During  the  steady  flight  condition,  the  airplane  is 

flying  with  wings  level  and  with  all  components  of  velocity 

zero  except  U and  W 
o o 

10.  Body  axis  system  is  assumed 


Dividing  the  force  equations  by  the  mass  m and  the 
moment  equations  by  the  appropriate  moments  of  inertia 
yields  terms  of  the  form: 

1 jX  u and  'i  AL  r.  < B1  ) 

m 0u  T ~Jr 

XX 

Replacing  the  first  term  by  X and  the  second  term  by  L, 

u r 

simplifies  the  notation.  Tnese  quantities  are  called 

"dimensional  stability  derivatives".  Using  this  short  hand 

notation  the  equations  of  motion  applicable  to  the 

longitudinal  problem  can  be  written  as  follows;  complete 

set  of  equations  is  presented  in  Hef.5) 


u + H q + qG  cos  0 =Xu*Xu+Xq+Xq 
o o u u q g 

♦ X w + X . w ♦ Xd  C?  + X<>  0 + Xy  o 

w w ee  ee  ee 


♦ Xd  Q * 6 * X£  a 

F F F F d B B B 


♦ *a  <? 

B 


♦ cos 


£ T u + cos  £ Tx  ($ 

u i °P.PM  RPM 


( B2  ) 
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w - (J  g ♦ g0  sin  0 = Z u ♦ Z.u  ♦ Z g * z.q 

o o u u g g 

+ Z * + Z.w  ♦ ZA  Q + Zi  <$  ♦ Z^  <5 

v v °ee  ee  e e 


♦ 


♦ 


6 

B 


■f 


6* 

B 


(sin 


(sin£  ) 


<5 

RPM 


( B3  ) 


q = a u + MtU  + M q ♦ flj  + H » + B,W  ♦ Q 

u u g g w w ee 


♦ 


♦ 


♦ 


+ 


♦ 


<? 

B 


* 


2 m 

MV  3 * — T u 

0 3 B I u 

yy 


z ro 

+ 1 T <y 

I ° EPS  RPM 

yy 


( 84  ) 


In  order  to  simplify  these  equations  it  was  assumed  that 
all  disturbances  were  small  such  that  the  higher  order  terms 
were  negligible.  Also  it  was  assumed  that  the  flaps  and 
speedbrakes  were  retracted  and  that  thrust  effects  were 
negligible.  Kith  these  further  assumptions  the  w eguation 
became: 

w - 0 g ♦ g6  sin©  = Z\i+Zq  + Zw  + Z*&  { B5  ) 

o o u g v we  e 

It  was  further  assumed  that  the  flow  was 

quasisteady.  Because  of  this  all  derivatives  with  respect  to 

the  rates  of  changes  of  velocities  were  omitted  with  the 

exception  of  those  with  respect  to  w which  were  retained  to 

account  for  the  effect  on  the  horizontal  tail  of  the 

downwash  from  the  wing.  If  a wing  moving  through  an  air  mass 

with  an  angle  of  attack oc  and  a steady  forward  velocity  U 

o o 
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is  suddenly  giver  a vertical  velocity  change  w without 
changing  the  forward  velocity,  the  angle  of  attack 
perturbation  is  as  shown  in  Figure  1. 


( Initiol) 


( Perturbation  In  w) 


Figure  1.  Angle  of  Attack  Due  to  Vertical  Velocity 

Perturb?  tion 


Thus  changes  in  angle  of  attack,  , are  proportional  to 
the  vertical  velocity  perturbation  w.  Thus  we  can  write. 


w ♦ W 

o 

tan  ( OC  ♦ oc  ) = * 

o VU2  ♦ <v  ♦ W )2 

o o 


< B6  ) 


anc  by  the  soall  angle  assumption 


and  CC 


( B7  ) 


Taking  the  derivative  yields. 


w = U a 
o 


( B8  ) 


also  it  can  he  ‘.  h vn  that. 


— = Liait 
w w — » 0 


and  substituting  equation  B 8 yields. 


iZ  1 lZ 

= Limit — — - 

aw  oc — ►o  U a 
o 


U 


( BIO  ) 


also  if  we  redefine  the  Z dimensional  derivatives  such  that, 
typically, 


* 1 

, z = — z 

q U q 


( Ell  ) 


we  can  now  substitute  into  equation  B 3 which  yields  the 
equation , 


CC  = z oc  ♦ 2 u + Z q 6 sin  0 + Za  (5 

cC  u q U o e e 

o 

i A 

vbere  z = 1 ♦ Z 

q q 


( B12  ) 


A similar  expression  was  obtained  for  u usinq  the 
simplifying  assumptions  and  the  expressions. 


i£ 


i 

Limit  

cC— 9 o U 

0 


aX 


1 

X 

U OC 
o 


and  the  expression  for  u was  therefore  written. 


u * X CC  ♦ X u ♦ X#  q -ge  cose  ♦ Ia  O 
CC  u q o e e 


( 313  ) 


( B14  ) 
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( B 15  ) 


where  X'  * (X  - oc  0 ) 
q g o o 

Finally  the  expression  for  q followed  directly  as, 

q = « 0 C u ♦ H g 6 

u gee 

The  model  equations  for  the  longitudinal  case,  including 
an  angle  of  attack  due  to  gusts  term,  were  written  in  oatrix 
fora  as. 


x « Fx  ♦ Gu  ♦ |"% 


( B16  ) 


where 
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x 


OC 

u 

9 

6 

OC 

9 


(angle  of  attack  perturbation) 
(longitudinal  velocity  perturbation) 
(pitch  rate  perturbation) 

(pitch  attitude  perturbation) 

(angle  of  attack  due  to  gusts) 
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Z' 
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Note:  Z ' in  the  aodel  is  equivalent  to  the  equation  of 

9 

■otion  ter  a (1  ♦ Z ) 
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6 

The  measurement  equations  written  in  matrix  form  were, 
y = Hx  ♦ Du  ♦ v ( B 1 7 ) 


1*.  K<X 

K ( *Ct  ) 4 q 
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Note;  Selected  stability  derivatives  are  primed  since 
they  differ  from  the  notation  used  in  Ref.  6,  however  they 
are  exactly  equivalent  to  the  derivatives  used  in  SCIDNT. 
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AIRPLANE  AND  GROUND  DATA  SYSTEMS 


A.  AIRPLANE  DATA  SYSTEM 

1 . General 

The  demonstration  F-14A  ship  number  6 BU.  NO.  157937 
Mas  the  test  airplatve  utilized.  A detailed  description  of 
the  airplane  instrumentation  system  is  contained  in  Refs.  7 
and  9.  The  airplane  was  modified  to  be  compatible  with  the 
Naval  Air  Test  Center  Real-Time  Telemetry  Processing  System 
(RIPS).  The  instrumentation  systems  and  measurements 
selected  for  use  in  the  test  airplane  were  primarily 
designed  to  satisfy  the  Flight  Test  requirements  fcr  the 
F-14A  Performance  Evaluation  and  Demonstration.  The  test 
airplane  was  not  entirely  a production  configured.  F-*14A  in 
that  some  miscellaneous  cockpit  components  not  required  for 
flight  testing  were  removed.  The  production  nose  radome  was 
replaced  by  a fiberglass  and  metal  radome  fitted  for  an 
instrumentation  noseboom  assembly.  Additionally,  some 

production  equipment  was  reworked  in  order  to  accomodate 
various  instrumentation  subassemblies. 

2.  No^e  jJoom  Instrumentation 

The  nose  boom  assembly  included  Angle  of  Attack  and 
Angle  of  Side  Slip  Vanes  and  a Pitot  / Static  Head. 
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3.  Cockpit  Instrumentation 


The  Forward  Cockpit  instrumentation  was  primarily 
pilot  oriented  and  included;  Angle  of  Attack  and  Sideslip 
Indicators,  Airspeed,  Altimeter,  Machmeter,  Center  of 
Gravity  Normal  Acceleration  and  an  Instrumented  Control 
Stick.  The  Aft  Cockpit  Instrumentation  included  left  and 
right  engine  BPil  (N1),  wain  Fuel  Flow,  A/B  Puel 
Plow,  Outside  Air  Temperature  and  the  Instrumentation  Master 
Control  Console. 

4.  Compartment  Instrumentation 

The  instrumentation  subcomponents  were  mounted 
throughout  the  airplane  in  various  available  locations 
including,  the  Nose,  Forward  Compartments,  and  Equipment 
Compartments.  The  bulk  of  the  instrumentation  equipment  was 
located  in  the  forward  equipment  compartments. 

5 . System  Description 

The  instrumentation  system  consisted  of  the 
following  subsystems; 

o PCM  Signal  Conditioning 
o PCM 
o Timing 

o Magnetic  Tape  Recording 
o Telemetry 

o Power  and  Distribution 
o Ground  Support  Eguipment 

The  system  had  the  capacity  to  record  on  magnetic 
tape  and  to  telemeter  to  the  ground-data-re Juction-ccmpiex 
all  the  PCM  data  for  real-time  analysis.  The  PCM  Signal 
Conditioning  Subsystem  consisted  of  active  pre-sample  filter 
chassis  and  some  special  conditioning  required  for  certain 
measurements.  The  pre-sample  filters  had  cut-off  frequencies 
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of  5,  10,  20,  40,  and  80  Hz,  The  PCM  subsystem  consisted  of 
an  analog  multiplexer,  analog  to  digital  converter,  digital 
multiplexer,  and  the  required  control  logic.  The  subsystem 
had  the  capacity  of  accepting  up  to  332  analog  parameters  at 
frequency  response  from  dc  to  40  Hz.  Discrete  data  (On/Off 
signals)  and  a limited  amount  of  digital  data  and  time  words 
were  time-di»'ision  multiplexed  into  the  output  data 
stream.  The  output  data  was  presented  as  a succession  of 
twelve-bit  (including  parity)  parallel  data  words  in  NRZM 
format  which  were  recorded  on  the  digital  tape.  A serial 
output  in  NRZL  format  for  entry  into  the  "L"  Band  telemetry 
subsystem  was  provided.  The  PCM  format  is  contained  in  Ref. 
10.  An  IRIG-t  Time  Code  Generator  was  installed  to  provide  a 
common  time  correlation  for  magnetic  tape  and  telemetry 
data.  The  Magnetic  Tape  Recording  Subsystem  consisted  of  a 
PCM  Recorder  which  had  the  capability  to  record  the  PCM  data 
and  Time  Code  information.  The  Digital  Recorder  contained  a 
sixteen  track,  IBIG  compatible  head  capable  of  storing  one 
hour  of  digital  data  on  one  inch  tape  at  a tape  speed  of  30 
IPS.  The  data  was  recorded  in  parallel  using  saturation 
recording  technigues.  This  portion  was  capable  of  bit 
densities  of  1000  bits  per  inch.  The  Telemetry  Subsystem 
consisted  of  a signal  combining  unit,  an  "L"  Band 
Transmitter  and  two  blade  antenna.  The  operating  frequencies 
could  be  tuned  to  coinside  with  those  in  use  at  various  test 
facilities.  The  entire  Flight  Test  Instrumentation  System 
was  powered  from  the  aircraft  primary  (generating) 
system.  Should  either  the  left  or  right  engine  or  generator 
fail,  the  instrumentation  system  would  continue  until  such 
time  as  both  generators  failed.  The  distribution  system 
consisted  of  a 28V  dc  transf ormer/rectif ier  (T/R/,  and  a 
power/control  distribution  box,  which  provided  electrical 
system  protection  and  isolation  from  the  aircraft  power 
systems  through  individual  circuit  breakers.  The  equipment 
required  to  perform  checks  on  the  installed  instrumentation 
systems  was  configured  in  mobile  test  carts.  These  carts 
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required  only  60  cycle  power  for  operation  and  provided  the 
operator  with  a tool  for  conducting  pre-flight  and  post 
flight  tests  on  the  tape-  recording  system. 


B.  GHOUND  DATA  SYSTEM 

The  Naval  Air  Test  Center,  Patuxent  River,  Maryland, 
employs  a Real-Time  Telemetry  Processing  System  (BTPS)  as 
universal  support  for  testing  on  several  types  of  aircraft 
and  weapons  systems, Ref.  10.  The  computer  system,  designed 
by  Xerox,  consists  of  a Sigma  9 general  purpose  computer  as 
the  central  processor  supported  by  three  Sigma  3 computers 
performing  the  data  input  and  output  processing.  RTPS 
provides  the  capability  to  process  two  complete  and 
independent  tests  (data  "streamsM)  simultaneously  from 
either  tape  or  telemetry  sources.  These  "streams"  can 
support  totally  unrelated  and  asynchronous  test  flights  or 
can  support  setup  for  one  flight  while  processing  another 
flight. 

Each  Data  Channel  accepts  up  to  512  channels  of  PCM  data 
and  32  parameters  of  FM  data.  This  data  is  processed 
(individually  or  mixed)  by  a Sigma  3 and  preprocessing 
Programmed  Algorithm  Unit  (PAU) . Up  to  50,000  values  of 
Engineering  Units  (EU)  data  per  second  can  be  fed  to  the 
Sigma  9 while  simultaneously  being  stored  on  digital  tapes 
for  later  playback.  Two  large  screen  interactive  graphic 
display  stations,  that  provide  real-time  "hands  on"  data 
analysis  and  presentation,  are  controlled  by  the  Sigma  3 
Display  Processor.  Batch  jobs  can  be  processed  when  the 
real-time  and  non  real-time  operational  modes  are  not 
utilizing  the  resources  of  the  system  to  100  percent.  The 
batch  job  is  defered  if  additional  demands  are  made  by 
either  non  real-time  or  real-time  tasks  during  batch  job 
processing. 

The  specific  flight  test  data  utilized  by  the  parameter 
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identification  program  were,  V^,  ALPHA  FRL,  pitch  rate, 

pitch  attitude,  NZ  NET,  and  average  horizontal  stabilizer 

position.  The  raw  data  (digital)  was  used  to  generate  an  ED 
tape  for  input  to  the  parameter  identification  program,  in 
generating  the  EU  tape  corrections  were  made  to  the  measured 
values  as  required  in  accordance  with  Ref  11.  For  example 
corrections  to  angle  of  attack  included  nose  boom  bending 
moments,  fuselage  bending  moments  and  pitch. 
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Manufacturing  labor  costs  refer  to  all  direct  labor 
necessary  to  machine,  fabricate,  and  assemble  the  major 
aircraft  structure.  This  includes  the  direct  labor  portion 
of  components  that  are  built  off-site.  The  labor  necessary 
to  install  these  components  and  purchased  components  is 
also  included.  Man-hours  required  to  fabricate  the 
purchased  parts  and  materials  are  excluded  (Ref  4:27-28). 

Manufacturing  materials  include  raw  and  semi fabricated 
materials  plus  purchased  parts  and  equipment  used  in  the 
construction  of  the  airframe.  Purchased  parts  are  general 
hardware  items  such  as  electrical  fittings,  valves,  and 
hydraulic  fixtures.  Purchased  equipment  is  other  general 
purpose  items  such  as  actuators,  motors,  generators, 
landing  gear,  instruments,  and  hydraulic  pumps.  Such 
equipment  may  be  purchased  by  the  contractor  or  provided 
by  the  government.  However,  if  a piece  of  equipment  is 
specifically  designed  for  a particular  aircraft,  it  is 
considered  subcontracted  and  excluded  from  manufacturing 
materials  (Ref  4:31). 
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